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Abstract 

We analyse the naturalness of the Minimal Supersymmetric Standard Model (MSSM) in 
the light of recent LHC results from the ATLAS and CMS experiments. We study non- 
universal boundary conditions for the scalar and the gaugino sector, with fixed relations 
between some of the soft breaking parameters, and find a significant reduction of fine- 
tuning for non-universal gaugino masses. For a Higgs mass of about 125 GeV, as observed 
recently, we find parameter regions with a fine-tuning of 0(10), taking into account 
experimental and theoretical uncertainties. These regions also survive after comparison 
with simplified model searches in ATLAS and CMS. For a fine-tuning less than 20 the 
lightest neutralino is expected to be lighter than about 400 GeV and the lighter stop can 
be as heavy as 3.5 TeV. On the other hand, the gluino mass is required to be above 1.5 
TeV. For non-universal gaugino masses, we discuss which fixed GUT scale ratios can lead 
to a reduced fine-tuning and find that the recent Higgs results have a strong impact on 
which ratio is favoured. We also discuss the naturalness of GUT scale Yukawa relations, 
comparing the non-universal MSSM with the CMSSM. 
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1 Introduction 



Recently, the LHC experiments ATLAS and CMS have reported the discovery of a new 
resonance, mainly based on excesses of events in the 77 and ZZ* channels. The combined 
significance amounts about 5a for both, ATLAS and CMS 1 1 . The new particle is compatible 
with the Standard Model (SM) Higgs boson with a mass of rrih ^ 125 -i- 126 GeV |2|. In the 
context of low-energy Supersymmetry (SUSY), such a value still lies in the predicted range 
of the Minimal Supersymmetric Standard Model (MSSM), provided quite heavy top squarks 
(with masses > 0(1) TeV) and/or large left-right stop mixing (see for instance [s]) are 
present. Regarding direct searches for SUSY particles |4|, based on ^ 5 fb~^ of data, only 
negative results have been reported. 

From the point of view of the MSSM, the negative results of direct SUSY searches and 
the comparatively heavy Higgs mass are well consistent with each other. However, as the 
main motivation for SUSY is the solution to the hierarchy problem, which naively requires 
SUSY masses of the order of the electroweak symmetry breaking scale, one may ask whether 
the MSSM can still accomplish this task in a natural way (or at least for a moderate fine- 
tuning price). Regarding specific MSSM scenarios, it turns out that in the constrained MSSM 
(CMSSM) a significant amount of fine-tuning (at least > (9(100)) is unavoidable to explain a 
Higgs mass of rrih ^ 125 -i- 126 GeV. Following our paper |5 on the CMSSM, we thus find it 
interesting to re-address the naturalness of the MSSM in the context of models where certain 
universality assumptions on the SUSY parameters are relaxed (e.g. for scalar and gaugino 
masses) . 

In order to envisage the most promising ways to reduce the fine-tuning, we are going to 
adopt the following strategy: we start with considering a generic setup with 17 independent 
parameters inspired by the phenomenological MSSM (pMSSM), in order to identify those 
giving the dominant contributions to the fine-tuning. This allows us to identify the rigid 
relations among the parameters which can decrease fine-tuning. After discussing several 
possibilities, we will mainly focus on certain non-universal relations among gaugino masses, 
whose possible impact on fine-tuning have been studied before the LHC results by |6,7|J^ In 
this context, we will study the interplay between the fine-tuning, the model predictions for 
the light Higgs mass and the Grand Unified Theory (GUT) scale ratios of the the gaugino 
masses and of the third family Yukawa couplings, which are key quantities for discriminating 
among different SUSY (and SUSY GUT) models. In our analysis, we include various relevant 
experimental constraints, e.g. from BR(6 57), BK{Bs ^^(Bu tVj-). We 

do not impose strict constraints from requiring that the {g — 2)^ deviation from the SM 
(currently at the level of 3.2(j [IT]) is explained by SUSY. 

The rest of the paper is organized as follows: in the next section we study in a semi- 
analytic way the dependence of the electroweak scale on GUT scale boundary conditions 
whose choice is inspired by the pMSSM. Thanks to this analysis we can identify regions with 
low fine-tuning for certain limiting cases and identify as most promising case non-universal 
gaugino masses. In section 3 we revisit briefly the dependence of GUT scale Yukawa coupling 
ratios on low energy supersymmetric threshold corrections. Thereafter we make an extensive 
numerical analysis of the fine-tuning in an MSSM scenario with non-universal gaugino masses 
at the GUT scale and compare the results to various experimental results, most importantly 

recent discussion of the fine-tuning price of a 125 GeV Higgs after the first hints in December 2011 
within several SUSY models has been also given in Is] (see also [gllTo]). 
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the recent discovery of a new resonance around 125 GeV compatible with a Higgs boson. In 
the final section 5 we summarize and conclude. 



2 Fine- Tuning in the MSSM 



Fine-tuning in the MSSM has been extensively discussed in the literature starting with [12] 
As we are going to follow the same approach as in our recent paper |5|, we give here just a 
brief summary of the most important points. 

In the MSSM, by minimizing the scalar potential, the Z-boson mass can be computed in 
terms of the supersymmetric Higgs parameter and the soft SUSY breaking mass terms of 
the up- and down-type Higgs doublets. For moderate and large tan/3, one finds at tree level 



(1) 



The value of the low energy observable Mz can thus be obtained in terms of the funda- 
mental parameters of the high-energy theory, by considering the renormalization group (RG) 
evolution of ii and m/f^, which are assumed to be both of the order of the SUSY breaking 
mass scale, from the GUT scale down to low energy. Since experimental constraints force the 
SUSY scale to be somehow larger than M^, a certain amount of tuning is needed. In order 
to quantify fine-tuning, the following measure has been introduced [12j 



d log Mz 




a dMl 


5 log a 




2M| da 



(2) 



reflects the dependence of Mz on the variation of a given GUT-scale Lagrangian parameter 
a. This definition not only encompasses the obvious fine-tuning in fi which is needed to fulfill 
Eq. Q, but it also covers the tuning needed to have a small rn^Hu^ while other soft SUSY 
breaking parameters entering its RG evolution are relatively large. The overall measure of 
fine-tuning for a given parameter point is then defined as the maximum of all the single A^'s: 



A — max Aa . 



(3) 



Turning to the contributions A^ of the single parameters, let us first notice that the fine- 
tuning in fi is rather special as its low-energy value is determined by imposing the correct 
Mz. Hence, setting aside RG effects for the moment, A^ can be approximately expressed in 
terms of m|r as 



A, 



m 



+ 1 



(4) 



In the numerical analysis, the running of fi is of course included. It turns out that in the 
constrained MSSM A^ is often the dominant term in A followed by the tuning in Aq, mo and 
Ml/2- Performing a scan of the CMSSM parameter space, we found in fs] that the fine-tuning 
in II dominates in 49(74)% of the cases, in 29(0)%, mo in 11(0)% and M1/2 in 11(26)% 
for the {g — 2)^ constraint applied at 5(2)cr. 

The discussion of fine-tuning for the other parameters is more involved as they enter Eq. 
([1]) only via the RG evolution of m|^^^. It is, however, possible to approximately express 
m|^^ as a polynomial in terms of the high-energy parameters of the theory. 



''For an extensive list of references, see for instance [8| |l3| . 
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In the following we will give approximate formulas for the dependence of Toj^^ on the 
assumed fundamental GUT-scale parameters at low energies, which can tell us for which 
choice of parameters we can expect a reduced fine-tuning. 

Before we come to this we like to note that fine-tuning measures should be taken with some 
caution. It is a matter of individual taste how much fine-tuning one accepts as "natural" and 
it also depends on the fine-tuning definition used. Thus we want to stress that our aim here 
is mainly to compare different parameter regions and only the relative fine-tuning difference 
between them would render a region more attractive to us. 



2.1 Fine- Tuning a la pMSSM 



To express rn^^ in terms of high-energy parameters, our choice of the independent parameters 
is inspired by the so-called phenomenological MSSM (pMSSM) [14J (but with the parameters 
defined at the GUT scale). In particular, we consider different scalar masses for the chiral 
superfields of the MSSM. The first two generation sfermions are assumed to be degenerate, 
while we allow a splitting for the third generation. The gaugino masses are taken to be non- 
universal. Concerning the trilinears we introduce three independent parameters: one for the 
up squarks, one for the down squarks and one for the charged sleptons. The main difference 
to the usual pMSSM is in the Higgs sector. Namely, instead of low scale /x and the CP-odd 
Higgs mass m^Q, we take the GUT scale value of rn^^ and m|^^ to be the free parameters 
and denote them by m^^ and m^^ respectively. 

The dependence of tti^^ on the GUT scale parameters can be deduced from its RGE (see 
e.g. [I5i) and can be written as 



(5) 



where the Ni = (Mi, M2, M3, Af^Ai^^A^) are assumed to be real and the matrix bij is a general 
symmetric 6x6 matrix. 

To make the deviation from the universality assumption of the CMSSM more evident we 
introduce the following dimensionless quantities: 



rj^ = Ma/Ms with a = 1, 2, 3 



rji Ai/Ms with i^t.h, r. (6) 

1 and we can 



Note that the convention we use here and in the rest of the paper implies 773 = 
recast the Ni as 

Ni = Ms - (771, 772, ^3, Vt, Vb, Vr) . 
The fine-tuning measure introduced in Eq. ([s]) can now be easily expressed as 

.2 



A. 



max A 2 , 



for the scalar masses, whereas for the other soft terms we find 



A 



(no sum over i), — max A 



Ni 



and 



A, 



Ml 



(7) 



(8) 



(9) 



(10) 
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After this general introduction to fine-tuning a la pMSSM, we turn to expressing m|^^ at 
the SUSY scale for a point which corresponds to M3 = 0.6 TeV, rui = 1.5 TeV for every z, 
r]i = Ni/Ms = (1, 1, 1, — 5, — 5, — 5), tan/3 = 30 and positive. This point is in agreement 
with the Higgs mass range quoted and the BK{Bs bound |16| and we find 

mjj^{MsusY) = -0.0459m|^ + 0.0988m|^ - 0.0469m|^ + 0.0488m|^ - 0.0541m|^ 

- 0.3347m^ - 0.2500m^~ - 0.0154m2^ + 0.0245m? - 0.0236m| 

W3 (^3 ^3 ^3 

+ 0.6481m|^ + 0.0273m^^ 

- M|(1.2865 - 0.0216r?i - 0.0242r?? + 0.0230r?2 - Q.2mril + 0.0813i?i%) 
+ M|(0.2521r?4 + 0.0208% + 0.0175r?^) 

+ M| 771 (0.0087r/t - 0.0028% + 0.0053r?r) 
+ M| r;2(0.0852i7t - 0.0086% - 0.0114r?^) 
+ M|(0.0022r7^ - Q.UUrj^ - O.OOOlr?^) 

+ M|(0.0068%% + 0.0017%r/^ + 0.0007%r?^) , (11) 

where r/s is set to be 1. One can already spot here the most important contributions to 
come from the stops, the gluinos and the Higgs sector itself. The coefficients aj appearing in 
Eq. ([5| can be directly read off from Eq. (11). In order to show the correlations among the 
parameters % 
matrix: 



Ni/M^ in a transparent way, we display the coefficients bij as a symmetric 



/0.0242 



-0.0813 
0.2177 



0.0216 
-0.0230 
-1.2865 



0.0087 
0.0852 
0.2521 
-0.1244 



-0.0028 
-0.0086 
0.0208 
0.0068 
0.0022 



V 



0.0053 \ 
-0.0114 
0.0175 
0.0007 
0.0017 
-0.0001/ 



(12) 



Even though the coefficients and hij were numerically obtained in a specific point of the 
parameter space, we checked that they provide a reasonably accurate estimate of Tin?Hu wide 
regions of parameter space. The corresponding uncertainty of the coefficients is estimated to 
be of the order of 10-20%. This is good enough for the qualitative discussion we present in 
the following subsections on possible strategies to reduce the fine-tuning. Nevertheless, the 
results we present in section [4] are based on a full numerical analysis, with the RGE evolution 
of all parameters computed in each point of the parameter space. 



2.2 Our Strategy 

In the attempt to find SUSY models with reduced fine-tuning we assume that the under- 
lying model of SUSY breaking predicts certain fixed relations between the SUSY breaking 
parameters at high energy (e.g. at the GUT scale). As a consequence, cancellations among 
different contributions in Eq. ( pTj ) can occur which lead to a reduced fine-tuning. We note that 
such a behaviour is a well-known property of models with universal scalar masses (like the 
CMSSM), for which the coefficients ai in Eq. (11) cancel almost completely]^ However, in the 



This property is known under the name of "focus point" 17 , for a recent quantitative discussion see, 
B.g. [5]. 
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light of the recent Higgs results the CMSSM unavoidably has a certain amount of fine-tuning 
either from quite heavy squarks or from a large A-teim (in maximal mixing scenarios). In 
this paper, we therefore consider the possibility that the underlying theory predicts different, 
non-universal (but fixed) boundary conditions at the GUT scale, namely non-universal scalar 
masses (NUSM) and non-universal gaugino masses (NUGM). Although such relations may 
not hold exactly in realistic models, they may guide towards more natural SUSY scenarios. 



2.3 Fine- Tuning from the Scalar Sector 

We now discuss possible NUSM scenarios. To start with, let us note that the situation of 
universal sfermion masses in the CMSSM features a significant reduction of fine-tuning, since 



an almost-complete cancellation occurs automatically for the largest contributions in Eq. ( 11 ), 

i.e. between the term with m? and the terms with m\ and m% . From the point of view of 

Qs U3 

naturalness, such a fixed relation between m? and ni?~ , m% is desirable. 

'^^ Q3 U3 

Let us now consider non-universal fixed ratios for the soft scalar masses motivated by 
unified theories: if we assume that at the GUT scale the gauge interactions unify to one 
single interaction then the soft masses for the various fields are not independent anymore 
because they are (partially) unified in common representations. Two very prominent paths to 
unification are SU{5) \i8j and Pati-Salam |19 , which can both be embedded in 50(10) |20j. 
For simplicity we assume that possible higher order GUT symmetry breaking corrections, 
which might induce splittings within one representation are negligibly small and we remind 
that we have set the first two generations to have the same soft SUSY breaking masses. 

• We begin with the case of SU{5), where we find for Toj^^ 

^k(^suSY) = -0.0012m|^ + 0.0019m|^ - 0.6083m|^ + 0.0091m|^ 

+ OMSlml^ + 0.0273m|^ (13) 
+ gaugino masses and trilinear terms , 

where rn^ stands for the masses of the tenplet of SU{5)^ which contains the coloured 
doublet, the up-type coloured singlet and charged leptonic singlet (m^ = = = 
m\ ), and im?~ for the fiveplet of 5^7(5), which contains the down-type coloured singlet 
and the leptonic doublet (m|^ = m| = ^|,)- We allow an arbitrary splitting from 
the Higgs fields from the other scalar mass parameters and also among each other. One 
can clearly see that in order not to strongly increase the fine-tuning one has to impose 
a fixed relation where m?- ^ m? holds, otherwise one has to consider the fine-tuning 
from each parameter separately, cf. Eq. ([s]). 

• We come now to the Pati-Salam case for which we find 

^k(^suSY) = 0.0029m? - 0.0022m?^ - 0.3102m| - 0.2890m?3 + 0.6754m^ 
+ gaugino masses and trilinear terms , 

where in this case / denotes the left-handed doublets {im?~ = m? = m~), f denotes the 

right-handed doublets {mf = mj- = m?~ = m? ) and h denotes the Higgs bi-doublet 

Ui Ui h/\ * 

(m^^ — — '^D- "^^^ conclusions are similar to the SU{^) case. The fixed relation 
one now has to impose is m~ mp^ m\. 
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• Finally, let us discuss the situation in 50(10) GUTs. Here one considers the soft mass 
terms for the matter fields m^g^ and ^igg and for the Higgs fields m^^ = m^^ = m^g 
(with possible m^^ ^ m^^ from D-term splitting), obtaining 

ml^jMsusY) = -0.0031m?6^ - {).bm2mlQ^ + 0.6754m2o ^^^^ 
+ gaugino masses, trilinear terms and D-term splitting terms , 

which implies that again, a fixed relation beyond the ones from the GUT itself would 
have to be imposed (between and ^igg) in order to strongly decrease fine-tuning. 

A few additional comments are in order. To start with, one can see that the soft terms for 
the first two families enter with small coefficients into m|r . This illustrates the well-known 
fact that these soft masses can be significantly larger than the third family ones, without 
paying a large fine-tuning price. Furthermore, in the above discussion we have focused on fixed 
relations from GUT structures. As we have already noted, with these fixed relations alone one 
would increase fine-tuning compared to the CMSSM case. From the above equations one can 
imagine additional fixed relations on top of these structures which could in principle reduce 
fine-tuning by effectively leading to a cancellation in the contributions to However, we 

will not investigate such cases in more detail here. 

Another relevant remark is that we have only discussed here "direct" effects on the fine- 
tuning, i.e. the effects on the fine-tuning measure for a fixed SUSY parameter point. We 
note that there are also "indirect" effects expected, i.e. when introducing a non-universality 
allows to avoid certain constraints on the SUSY parameter space and makes regions with 
lower fine-tuning accessible (e.g. a "compressed" spectrum helps to evade LHC constraints). 

Finally, let us notice that small splittings in models where the soft parameters are universal 
(CMSSM-like) at the leading order and get only corrected by small additional contributions 
(e.g. by higher dimensional operators in flavour models, see, for instance, |21 ) can have a 
certain degree of non-universality without increasing the fine-tuning too much. However, 
large non-universalities (0(1)) again lead back to the same situation as in the general case a 
la pMSSM. 

In summary, we conclude that it seems unavoidable to have the scalar mass parameters 
- at least the high energy stop and Higgs soft masses - almost degenerate to keep the fine- 
tuning small, like for instance in the CMSSM. Non-universal scalar masses tend to make the 
fine-tuning worse rather than improving the situation. We will therefore not study this case 
numerically in more detail. 



2.4 Fine- Tuning from Gaugino Masses and Trilinear Couplings 

Now we turn to the gaugino masses and the trilinear SUSY breaking couplings. They appear 
together in the polynomial and hence should not be discussed independently. For the 
sake of simplicity we assume here that the trilinear couplings are universal rjt — rjjj — rj^ — rjA^ 

^ In principle, one can reduce the fine-tuning as well by assuming definite relations between 77^, 775, and rfr 
just as we will demonstrate here for the gaugino mass parameters. However, the coefficients involving 775 and 



?7r in Eq. ( 11 ) are very small and hence the values of 77^ and 77^ should be very large. 
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Figure 1: The two ellipses in the 771-772 plane which minimise the fine-tuning. The red ellipse 



corresponds to Eq. (18) and the black one to Eq. (19). 



In this case the rnj^^ polynomial reads 

^k(^suSY) = -M|(1.2865 - O.O2I67/1 - 0.02427/? -f 0.02307/2 - 0.21777/1 -f 0.08137/1 7/2) 
+ M|7/a(0.2904 + 0.01127/1 + 0.06527/2) - 0.1131M|7/i 
+ scalar masses 

= + /2(r?i, %)r?A + /3r?i)M| + . . . , (16) 

where we have again set 773 = 1. Now the question arises for which values of the parameters 
the fine-tuning, as defined in Eqs. ^ and (10), is minimised. 



The overall fine-tuning, considering only this sector is 
A 



Mi 



Ml 



"2/l(??l,??2) + /2(??1,'?2)'?a| ' 



Setting these two equations to zero we find two possible solutions: 

Solution 1: 172) = 0, ??A = , (18) 

Solution 2: /,(,„,2) - ,A - • (19) 



With the parameters in Eq. (16), the two solutions define two difiFerent ellipses in the 7/1-7/2 
plane, see Fig. [l] As we will see in our numerical analysis later on this is a generic feature of 
the parameter space. 

As the A-terms are crucial for the prediction of the physical Higgs mass, it is interesting 
to discuss the dependence on rfA in more detail. For the first solution this is trivial, while the 
second solution can be written in terms of 7/1 and 7/2 depending on tja so that we find 



7/1 = -8.0849 + 6.07577/A T 4.9510^(2.1354 - 7/A)(r/A - 0.0962) , 

7/2 = -3.0681 + 2.42887/A ± 0.8483^(2.1354 - 7/A)(r/A - 0.0962) , ^^^^ 

which gives real ratios only for 0.0962 < 7/^ < 2.1354. Interestingly, this implies for example 
that in our benchmark point with 7/^ = — 5 the fine-tuning cannot be made arbitrarily small 
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by choosing appropriate rji and r]2- Still the fine-tuning can be significantly reduced: for this 
benchmark point we find a minimal fine-tuning in the gaugino sector of 



A = 4.95 



Mi 



(21) 



for 771 -11.7 and 772 -4.6 compared to A = 7A9\M^/M'^\ for 771 = 772 = 1. 

Note that we have discussed here the fine-tuning in M3 and 77^. As we mentioned above, 
there are regions of the parameter space where the fine-tuning in is dominant. Nevertheless, 

is also reduced by the choices of 771, 772 and tja we have just discussed. From Eq. Q we 
see that making l^nf^J small reduces the fine-tuning in fi. For solution 1 this is obvious. In 
the case of solution 2, it is easy to check that the contribution to |^^^| from the gaugino 
masses and trilinears is set to zero, too. 

The first solution for the case 77^ = was already found in |7| - albeit with a slightly 
different notation setting 772 = 1 instead of 773 = 1 and thus giving hyperbolas instead of 
ellipses. The second solution, however, was not discussed therein. Instead, in |6| a related 
case was studied where also a fixed relation between the gaugino and the trilinear mass scales 
was assumed, while in our analysis 77^ is a free parameter. We will discuss which GUT 
scale gaugino mass ratios are favoured, comparing numerical with semi-analytical results, in 
section [H 



3 GUT Scale Yukawa Ratios and SUSY Threshold Correc- 
tions 

In GUTs quarks and leptons are unified in representations of the GUT symmetry group, which 
implies that the elements of the Yukawa matrices are generally related. As a consequence, 
fixed ratios of quark and lepton masses are typically predicted by GUT models. These ratios 
are given by group theoretical Clebsch-Gordan factors and hold at the GUT scale. For the 
third family, phenomenologically interesting possibilities are, e.g., the standard bottom-tau 
Yukawa unification yr/vb = 1 or the recently proposed ratio yr/Vb = 3/2 [22]. The GUT scale 
Yukawa ratios are thus interesting discriminators between GUT models. We are therefore 
also interested in which GUT-scale Yukawa coupling relations might be preferred due to 
lower fine-tuning. 

The phenomenological viability of GUT scale Yukawa ratios depends on SUSY threshold 
corrections |23|. Here we give some short comments on this important class of one-loop 
corrections to the Yukawa couplings of the down-type quarks and the charged leptons in the 
MSSM, which are sizeable if tan /3 is large. These corrections depend on the SUSY breaking 
parameters and hence the GUT scale Yukawa coupling ratios depend on the SUSY spectrum. 

For example for the bottom quark mass we can write 

rrib = VbVdi^ + ^b tan /3) , (22) 
the correction 65 being given by [24[j26] 

€6 « + €^ + €^ + e^ (23) 



8 



where 



2as II , . 



1 



6 n 



w 



167r2 
1 3^^ r?2M3 



H2{vq^,xi) + 2H2 {vj^ ,Xl] 



/2 



167r2 2 n 



H2{Vq^,X2) 



e" — 



H2{vQ^,Vu^] 



Here u j 



rnVMl vf 



function H2 reads 



WJ 



H2ix,y) 



xlnx 



+ 



infMi)//,^ for i 
ylny 



(24) 
(25) 

(26) 
(27) 

1, 2 and the loop 
(28) 



{l-x){x-y) {l-y){y-x) ' 

In the CMSSM one can neglect in a first approximation and because they are suppressed 
by the small gauge couplings. In NUGM scenarios this is in general not true anymore because 
the suppression might be compensated by an enhancement of the bino or wino mass parameter 
compared to the gluino one. In our example point we have seen that, for example, 771 = —11.7 
and 772 = —4.6 is preferred to have low fine-tuning. 

Without some knowledge of the parameter space or simplifying assumptions a quantitative 
statement is hence quite difficult. But in general one can expect significant corrections up to 
50 % for the GUT scale Yukawa coupling ratios (see, e.g. (251I26)). 



4 Numerical Analysis 

To improve our understanding based on the semi-analytical treatment above we now turn 
to our numerical analysis. To this end we employed a modified version of soft SUSY |27| to 
calculate the SUSY spectra, GUT scale Yukawa couplings and the fine-tuning. The modifi- 
cations take into account the implicit Mz dependence on the Higgs vev [28^ and the SUSY 
threshold corrections for the Yukawa couplings of all three fermion generations [23- - -26| . 
In order to compare with experimental constraints the observables BR(6 57), BK{Bs 



fi^ fi ) and BR(S^i rur) were calculated using Superlso 29 . The experimental ranges 
we considered for such observables are: BR(6 57) = (355 ± 24 ± 9)10"^ (SO), BR{Bs 
IJL-^IJL-) < 4.5 • 10-^ at 95% CL [16 and BR(S^ ^ rur) = 1.41 ± 0.43 |3l). Additionally, 
LEP bounds |32 on sparticle masses were applied and we discard points corresponding to a 
f LSP and color and charge breaking (CCB) vacua. We present results with and without the 
Higgs boson mass constraint included to study its impact. Applying constraints from direct 
LHC SUSY searches is not straightforward and we will only show a comparison to simplified 
models for illustration. 

In our opinion the 3.2cr level (^ — 2)^ discrepancy from the Standard Model expectation |TT] 
is an experimental evidence that, while certainly interesting, still requires a full experimental 
confirmation from next generation experiments, such as those proposed at J-PARC [33j and at 
Fermilab |34|. For this reason we have not imposed constraints from the anomalous magnetic 
moment of the muon. In fact in this scenario we can only find very few and isolated points 
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Figure 2: Lowest fine-tuning in the (a) r]i-r]2^ (b) m-Ma, (c) tja-Ms and (d) Vt/yb-yr/yb 
planes consistent witii the experimental bounds described in the text. For explanations of 
the parameters see main text. The Higgs mass constraint is not included. The la errors 



on the quark masses 32 are taken into account by scaling the data points in the last plot 
correspondingly. 



which would satisfy {g — 2)^ at 2a and predict a Higgs boson mass in agreement with the 
recent results, all of them corresponding to a fine-tuning larger than (9(100). 

Moreover, we have not imposed constraints from neutralino relic density nor from direct 
and indirect dark matter searches. Such constraints can be evaded assuming a non-standard 
cosmological history or a different LSP, like the gravitino or the axino. Nevertheless, it is 
interesting to note that for 772 ^ 0.5 771, the lightest neutralino is dominated by its Wino (or 
Higgsino) component, which implies that the relic density is strongly suppressed. In these 
regions, on the one hand, one cannot explain dark matter by thermal relic neutralinos, but, 
on the other hand, there is no overproduction of neutralino dark matter. 
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4.1 Before the latest LHC Higgs and SUSY Results 

As previously discussed, we expect to find an ellipse shaped region in the 771-772 plane where 
fine-tuning is significantly lower than in other parts of the parameter space, especially when 
compared to the CMSSM. To study this in detail, we scanned the NUGM parameter space in 
the following region: the soft scalar mass parameters were assumed to be universal rrii — m 
for every i and varied from to 4.5 TeV, the GUT scale gluino mass M3 from 0.15 to 2 TeV, 
and 77t = 775 = 77t- = 77^ from —20 to 20 for tan/3 = 2, 10, 15, 20, . . . , 60. For the ranges of 
the gaugino mass ratios we took —40 < 771 < 45 and — 10 < 772 < 10. Both signs for ji were 
allowed. In the scan we have dropped all points with a fine-tuning larger than 200. 

Fine-tuning was calculated for the parameters 771, /x, — tjaM^ and M3. Nominally 
existing tuning in the ratios 771, 772 was neglected as we assume such relations to be rigid and 
given by some underlying model. We note as a side remark that in all plots where contours 
are shown instead of only points, outliers that are very far away from their neighbours are 
still shown as isolated points instead of being incorporated into the contour. In addition, 
unless stated otherwise, all plots show values marginalised over parameters not shown. 

As we can see from Fig. [2^, we indeed find an ellipse that allows for quite low fine- 
tuning and evades the experimental bounds. We see that the semi-analytical results from 
Fig. [l]do not very well reproduce the numerical results but still they give a good qualitative 
understanding. Compared to the CMSSM with the same experimental bounds, we find a 
fine-tuning A lower than 3 instead of just below 10 |5|. In this section, we will focus on this 
region of low fine-tuning and only consider points with an (771,772) combination lying on the 
ellipse defined by Amin < 10,^ An approximate analytical formula for the ellipse in Fig. [2^ 



is given by 

(771/15.0)2 + (772/2.6)2 = 1 . (29) 

Based on this, we take a closer look at the fine-tuning dependence on the other parameters. 
In Fig. [2)3 we can clearly see that introducing non-universal gaugino masses can significantly 
weaken the dependence of Amm on the gaugino mass scale compared to the CMSSM. While 
in the latter a fine-tuning of 20 is only possible for M3 < 400 GeV now we can go up 
to M3 ^2.2 TeV. On the other hand, the fine-tuning dependence on m (resp. mo in the 
CMSSM) is very much the same as expected. 

Furthermore from Fig. ^ we can roughly see the same behaviour for large ija as in the 
CMSSM [5], while for small ija from -2 to 4 we find two peaks. They lie at 7/^ = and ija ^ 2, 
which are the values for r]A required for solution 1 and the upper end of the range quoted for 



solution 2, cf. Eqs. ([18) [19]). 

Finally, concerning the GUT scale Yukawa coupling ratios shown in Fig. [2]l we find no 
qualitative difference compared to the CMSSM, besides the rather obvious fact that the fine- 
tuning price of all ratios is greatly reduced. The comparison of the b-r Yukawa coupling ratio 
of 3/2 with the case of b-r Yukawa unification, however, is slightly in favor of 3/2 with a 
minimal fine-tuning of A ^ 5 vs. 9. 
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planes consistent with the experimental bounds described in the text. For explanations of the 
parameters see main text. In comparison to Fig. |2]we have included the CMS experimental 
constraint rrih = 125.3 ±0.6 GeV [l] and a theoretical uncertainty of ±3 Ge V [38| for the Higgs 
mass calculation at each data point. The la errors on the quark masses [32 are taken into 
account by scaling the data points in the last plot correspondingly. 
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Figure 4: Dependence of gaugino mass contribution functions /g, /a (left plot /black, 
left plot/red, right plot respectively) to left-handed stop soft squared mass m^, right-handed 



stop soft squared mass m| and stop trilinear coupling At as defined in Eq. 30 For the ellipse 
of low fine-tuning the rough form of Eq. (29) was used and we have defined (j) = arg(7^i -\-ir]2). 
The corresponding value of fx for = 7^2 = 1 is shown as dashed line in the same color. Note 
that fu has been rescaled by a factor 5 for illustration. 



4.2 Results including the latest LHC Higgs and SUSY Searches 

In the light of the discovery of a possible 125-1-126 GeV Higgs boson in current LHC searches, 
it is interesting to re-analyse the consequences of non- universal gaugino masses for naturalness 
of the (non-universal) MSSM. In Fig. [s] we show the same plots as in Fig. [2] with the additional 
experimental constraint from CMS, ruh = 125.3 ± 0.6 GeV |1|, on top of which we include a 
theoretical uncertainty of ±3 GeV |38] for the Higgs mass calculation at each data point. 

As we can see from Fig.|3^, non- universal gaugino masses can accommodate this additional 
constraint even with a fine-tuning lower than 20 (actually just above 10 and in a very small 
region of parameter space even slightly below 10). This happens e.g., around 772 ^ 0, 771 ^ 15 
and to a lesser degree also around 771 ^ 0, 772 ^ 2.6 and in the intermediate part of the ellipse. 
The reasons why these regions are favoured after including the Higgs results can readily be 
seen from the beta functions flS] of the stop soft masses masses and the stop trilinear coupling 
at the GUT scale: 

/32 2 \ 

IGir^Pm^, D -gluTMi i- + 6r?i + -Tjjj = -y^uT Ml fQ{m, m) , (30a) 

/32 2 \ 

IQn^P^^, D -gluTMi i — + —vl) = -9GUT Mi fu{vi,m) , (30b) 

/ 22 26 \ 

167t^13a, D -^gut^3 f y + 67/2 + Y^Vij = "^gut^s fA^m.m) . (30c) 

The functional behaviour of /g, fu and with 771 = cos(/), 772 = r^y sin(/) and by means of 
the approximate form of the numerically obtained ellipse of Eq. (29), is shown in Fig. [4j It 



can be easily seen that the gaugino mass contribution to right-handed stop masses is greatly 
enhanced for the angles = 0, tt, i.e. 772 = 0, |77i| ~ 15, while it is not significantly away from 
the CMSSM value elsewhere. The contribution to the left-handed stop mass on the other 



Depending on the other parameters these points can have a larger fine-tuning but, interestingly, this does 
not significantly change the appearance of all affected plots, indicating that being on the ellipse is indeed a 
necessary condition for small fine-tuning. 
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Figure 5: Lowest fine-tuning in the mh-yr/yb plane consistent with experimental bounds 
described in the text for non-universal gaugino masses (left) and universal gaugino masses 
(right). The horizontal lines correspond to the la error of 0.6 GeV around 125.3 GeV as 
claimed by the CMS collaboration 1 . The la errors on the quark masses |32 are taken into 
account by scaling the data points correspondingly. In addition, a theoretical uncertainty of 3 
GeV 38 for the Higgs mass calculation at each data point is taken into account. 



hand receives a smaller but still sizable enhancement over the CMSSM value, but does not 
vary much along the ellipse. Additionally for = 0, i.e. 771 ^ 15 and 772 = 0, the running 
contribution to At receives a sizeable enhancement, which is present to a lesser degree also for 
= 7r/2, i.e. 771 = and 772 ^ 2.6. For angles in the other quadrants, where at least one 77^ is 
negative, however, the cancellation with the large gluino mass running contribution prevents 
from a similar enhancement in /a- However, in practice, these three effects turn out to be 
effective only when they work in unison to obtain a high enough Higgs mass. 

We note that the changes in Fig. [8)3 are not very suprising: regions with low masses get cut 
off. Also the changes in Fig. [Sj^ are not unexpected, the solution with 77^ > is disfavoured 
more than the solution with 77^ ^ 0, as the latter does not suffer from the cancellation with 
the gluino mass contribution to the running of the top trilinear soft term. Unsurprisingly, 
it is shown that either large 77^ or large M3 (hence large top trilinear coupling) is needed to 
accommodate a large Higgs mass. 

Concerning the Yukawa coupling ratios shown in Fig. [3]l we can see that the minimal 
fine-tuning required to get to a unified b-r Yukawa coupling ratio suffers more from the 
requirement of a consistent Higgs mass than the ratio of 3/2 does. Namely after applying the 
cut the fine-tuning for yr/Vb = 3/2 can go down to A = 30, while yr = yb requires at least 
a A of 60. This is also nicely illustrated in Fig. [5] where we have shown the interplay among 
the Higgs mass, the GUT scale yr/yb ratio and the amount of fine-tuning. 

Another important aspect of current LHC experiments are the searches for supersymmetric 
particles. Since changes of the gaugino mass ratios 771 and 772 can significantly alter the 



composition of the lightest neutralinc^^ as well as the mass splittings controlling jet energies 



and missing Et from the cascade decays, the exclusion regions for the CMSSM do not apply 



Actually, low fine-tuning generally also means quite low /i, so Higgsino-like lightest neutralinos and 
charginos are likely. 
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Figure 6: Lowest fine-tuning in various planes used for simplified models. Only points 
consistent with experimental bounds described in the text are shown. In addition, for the plots 
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Table 1: Selected ratios and the minimal possible fine-tuning they allow after requiring the 
experimental constraint rrih = 125.3±0.6 GeV 1 and a theoretical uncertainty of ±3 GeV |38 
for the Higgs mass calculation. Only ratios that can reduce the fine-tuning by at least 50% 
compared to the unified (CMSSM) scenario are shown. For more details on the origin of these 
ratios, see, e.g. [7^ |37(|39| . The results are illustrated graphically in Fig. [sj 



anymore. A full event and detector simulation would, however, go beyond the scope of this 
study; thus we make use of exclusion bounds derived in so-called simplified models |4|[36|. For 
simplicity we just compare the spectra found by the numerical scan with the most stringent 
bounds in several kinds of simplified models. While this is not certainly a rigorous approach, 
it should give a good feeling how endangered by exclusion each point is. 

The resulting situation is shown in Fig. [6j As we can see in most cases only the region 
with A < 3 is partially inside the excluded region, while even A < 10 extends far beyond the 
bounds. Requiring a Higgs boson within la of the experimental measurement [T] (including 3 
GeV theoretical uncertainty as explained before) excludes even more than what direct SUSY 
searches do in some areas. It is thus important to note that, even in this quite restrictive 
approach, all of the parameter space with high ruh and A < 20 is safe from being excluded 
by the LHC in the near future. 

For future searches it is also interesting to look at the correlation between the Higgs and 
several sparticle masses, which is shown in Fig. [7[ There we can see that for low fine-tuning 
relatively light neutralinos and charginos are expected. This is due to the fact that /i is 
small in this region, so that we expect Higgsinos to be light. We also see that the gluinos 
and the lightest stop are rather heavy in our scenario with low fine-tuning (A < 20 for 
rUg ^1.5 TeV and m^^ ^1-0 TeV). The LHC still did not reach this parameter region and 
hence the statement that the natural MSSM parameter space is already ruled out seems to 
be premature. 



4.3 Favoured Non-Universal Gaugino Mass Ratios 

Fixed non- universal gaugino mass ratios may originate from various high energy models, for 
instance from GUTs or orbifold scenarios (see, e.g. [7,37,39] for discussions). Tab.[l]and Fig. 
[8] show examples of proposed fixed ratios 771,772 where we find that the fine-tuning can be 
reduced by more than 50% compared to the CMSSM. 

Interestingly, among the orbifold models G-H of (7,39 , from the full numerical results 
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Figure 8: Lowest fine-tuning in tiie r]i-r]2 plane consistent the experimental constraint ruh = 
125.3 ±0.6 GeV 1 and a theoretical uncertainty of ±3 GeV 38 for the Higgs mass calculation 
at each data point. Theoretically motivated ratios that reduce the fine-tuning compared to the 
unified scenario by at least 50% and the analytical expectation are also included. The CMSSM 
ratio is marked as empty diamond, circles correspond to ratios derived from GUT symmetry 
breaking [37^ and squares to ratios found in the so called O-II model in [7{[39]. For details see 
Tab.[l] 
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including the constraints from the latest Higgs results we find that the option (771,772) = 
|) from Sqs = —6 has the lowest possible fine-tuning (Amm = 13), whereas in jrj, 
based on analytical estimates before the Higgs results were available, the preferred ratio was 
(r?i,r/2) = (f,3)from5GS = -5g 

The GUT ratio with lowest fine tuning (Amin = 12) turned out to be (771, 772) = (10, 2). For 
the ratios found to be favoured in |7|, we find significantly higher fine-tuning, e.g.: Amin = 82 
for {m,V2) = (-5,3), A„i„ = 141 for (771,172) = (-^,-1), A„i„ = 143 for (171,772) = (1,-i). 

As it can be seen as well from comparison of Fig. [2^ and [3^, the inclusion of the Higgs 
results has strong effects on the favoured non-universal gaugino mass ratios. 

Finally, as discussed at the beginning of section [4| for 772 < 0.5 771 the neutralino is dom- 
inated by its Wino (or Higgsino) component which implies that the relic density is strongly 
suppressed and there is therefore no danger of overproducing thermal neutralino dark matter. 
Among the ratios listed in Tab.[l]this applies to (771,772) = (10,2) and (771,772) = (^, 1). For 
the other ratios of Tab. [TJ 772 > 0.5 771 holds and, at least in principle, thermal neutralino 
dark matter could be possible (for a standard thermal history of the universe). As already 
mentioned, we have not applied dark matter constraints in our analysis since they depend on 
the thermal history of the universe. 

5 Summary and Conclusions 

We have re-addressed the question of the naturalness of the MSSM, in the light of the discovery 
of a Higgs-like resonance recently announced by the LHC experiments |1|. We focused on 
models with non-universal boundary conditions of the SUSY breaking parameters at high 
energy and compared them with the CMSSM. Our basic assumption was that the departures 
from universality are a consequence of an underlying mechanism, possibly associated to SUSY 
breaking or GUT dynamics, giving fixed relations among the parameters. In order to identify 
which of such relations can lead to a reduced fine-tuning, we first considered a general high- 
energy parametrization of the soft terms. Then we discussed scenarios with non-universal 
scalar and gaugino masses. We found this latter possibility particularly promising (as already 
noticed, e.g., in (6|[7]) and we thus studied it in detail, computing numerically the fine-tuning 
measure A. 

We found that, considering the uncertainty related to the theoretical prediction, models 
with non-universal gaugino masses can account for a Higgs mass in the CMS la range 125.3 ± 
0.6 with a fine-tuning price of 0(10), in contrast to the CMSSM that requires Amin ^ O(IOO). 
Thus the MSSM with specific gaugino mass ratios still represents a comparatively natural 
scenario. In Fig. [2] and |3} we have shown the values of 771 — Mi /Ms and 772 — M2/M3 giving 
a low fine-tuning before and after applying a constraint on the Higgs mass. Interestingly, 
some of the ratios discussed in the literature lie in (or are close to) the new low fine-tuning 
region, while others are disfavoured when the new Higgs results are included. Including the 
Higgs results we found that particularly favoured ratios (with Amin = 0(10)) are now, e.g., 
(^1, V2) = (10, 2) which may originate from SU{5) GUTs and (771, 772) = (^5 i) from orbifold 
scenarios of type O-II with Sqs — —6. 

Furthermore, allowing for non- universal gaugino masses at the GUT scale, we have ana- 
lyzed the fine-tuning price of different values of the GUT-scale ratio yr/Vhi that represents an 
important handle to discriminate among different GUT models. For rrih ^125 GeV, we found 

^^fcs is a negative integer constant associated with Green-Schwarz anomaly cancellation (cf. [7|p9]). 
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that b-r Yukawa unification corresponds to A > 60, while the alternative ratio yr/Vb — 3/2 
can be realized at the price of A > 30 only. 

Concerning the SUSY spectrum that is favoured by naturalness and rrih ^ 125 GeV, 
we found that for the least tuned data points with fine-tuning A less than 20 the lightest 
neutralino is expected to be lighter than about 400 GeV and the lighter stop can be as heavy 
as 3.5 TeV. On the other hand, the gluino mass is required to be above 1.5 TeV. Comparing 
the predicted spectra with the LHC exclusions derived in a set of simplified models, we could 
conclude that the regions of lowest fine-tuning are at present only poorly constrained by direct 
SUSY searches at the LHC. 

Let us finally remark that, although the CMSSM is certainly challenged from the fine- 
tuning point of view by a Higgs mass of ruh ^ 125 GeV, more general realizations of the 
MSSM, like the examples studied here, can still provide a relatively natural solution of the 
hierarchy problem and will probably require several years of data taking before being fully 
tested at the LHC. 
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